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I. INTRODUCTION 



Since its realization in dilute atomic gases, superfluidity of alkali atoms has been studied 
extensively in various externally controllable environments [l]. As a controllable parameter, 
the external trapping potential of a cold atomic system plays an important constrained role 
in current experimental setups. Most of the experiments are carried out in either magnetic 
or optical harmonic traps. A negative, but small quartic term is always present with the 
Gaussian optical potentials in current experimental setups. For the case of rapidly rotating 
Fermi gasses in harmonic traps, the force due to the trapping frequency almost balances the 
centrifugal force and the atomic cloud spreads in the plane perpendicular to the rotation 
axis. Therefore, an added positive quartic potential term ensures the stability of the fast 
rotating regime {2, 5|. Adding a positive quartic term is very important for a two-component 
superfiuid Fermi system where the system is expected to be in the quantum Hall regime at 
the limit of extreme larger rotations Q]. 

Although the interaction between the atoms is extremely important in a system, the 
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problems are made tractable and the essential physics is retained by assuming fully polar- 
ized non-interacting atoms. Experimentally such a system can be prepared by trapping a 
single component (single hyperfine spin state) of atomic species. Because of the large spin 
relaxation time of the atoms compared to the other experimental time scales, single species 
can be retained over the entire time of the experiments. For a system of fully polarized 
Fermi gas, the s-wave scattering is prohibited due to the Pauli exclusion principle, while 
higher order partial wave scattering is negligible at low temperatures. Therefore, we model 
the system as a collection of non-interacting atoms obeying Fermi-Dirac statistics confined 
by quadratic and quartic potentials given as 



= br^ + b,z'^ + -fr^ (1) 

where M is the atom mass, Wj's are the harmonic trapping frequencies {i = r,z), and 
= + y^. Notice that we are using the cartesian coordinate system with coordinates 
{x,y,z). We use Thomas- Fermi semiclassical approximation to derive various low tempera- 
ture thermodynamic properties as a function of anharmonic confinement 7 for both positive 
and negative values. Fast growing progress in the current experiments to manipulate quan- 
tum many body states in rotating Fermi gases motivated us to study the effect of positive 7 
values [5I. Experiments in harmonic traps are limited by the loss of confinement. This loss is 
due to centrifugal force when the rotation frequency approaches the axial harmonic trapping 
frequency. One can overcome this loss by adding a quartic term to the trap potential as 
already implemented in experiments 5|]^. On the other hand, the small negative quartic 
term present with the Gaussian optical potential can be canceled by adding an extra quartic 
term. We will show later, that the minimum negative 7 value at which the atomic cloud 
is confined depends on the total number of atoms present in the system. Thermodynamic 



p.openi. of Fe.n,i ..e. — i„ a h..onic ..p W .ee„ ..d« . Ret. fl. In a 

recent work, interplay between rotation and anharmonicity of a Fermi gas has been studied 
extensively in Ref. [t]. In the superfluid phase of a two-component gas, the effect of anhar- 
monicities on breathing mode frequencies in the Bardeen-Cooper-Schriefer- BEG crossover 



^ We do not consider rotation in the present work. Fast rotating Fermi gas can be linked to one of the very 
interesting phenomena; integer quantum HaU effect. 



3 

regime has been presented in Ref. Q|. Thermodynamic properties of a Bose gas confined in 
a quartic potential trap can be found in ref. ^\. 

This paper is organized as follows. In section II, we present the derivation of thermo- 
dynamic quantities within our framework of semiclassical approximation. In section III, we 
present our results discussing how these thermodynamic quantities depend on trap potential 
parameters. Finally in section IV, we draw our conclusions. 



II. FORMALISM 



We consider a polarized Fermi atomic system trapped in an external potential which has 
both harmonic and radial quartic components given in Eq. [H For the atoms at thermal 
equilibrium, the average number of atoms in the single particle state \a > with energy is 
given 



where fi is the chemical potential and /3 = l/{kBT), with ks the Boltzmann constant and 
T the temperature. The average number = J2a of atoms of the system fixes the 
chemical potential. The total energy U = J2a Q^a^a, where the spin degenerate factor g = 1 
because the fermions are polarized. When the number of atoms in the system is large and 
the average potential energy of atoms in the trap is much smaller than the kinetic energy 
of the atoms, Thomas-Fermi semiclassical approximation can be used. The semiclassical 
approximation means that, instead of Cq, we can use the classical single-particle phase space 
energy e(r, p) = p^/(2M) + V{r,z). In general, at high enough temperatures, the particles 
are classical so that the semiclassical approximation is valid. In the present paper, we are 
considering a spin polarized Fermi system. As a result of the Pauli exclusion principle, 
when the temperature goes to zero, Fermi atoms fill every energy level up to the Fermi 
energy Ep (the energy of the highest occupied state). If the system has large number of 
atoms, then the Fermi energy Ep of the system is larger than the single atom ground state 
energy^. Therefore, a spin-polarized Fermi gas can be well described in the semiclassical 



^ For a system with large enough atom numbers, majority of atoms will have average energy in the order 
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approximation even at very low temperatures. This is drastically different for bosons where 
thermal energy /c^T has to be compared with the ground state energy. 

Using the quantum elementary volume of the single-particle phase space, /i^, where h is 
the Plank constant, the total number of atoms and the energy of the system have the forms 



^ g_ f d?rd^p 

h^J Z-ie/^^^'f) + 1 ^ ^ 

and 

TT 9_ f e{r,p)d^rd?p 
h^J Z-ie/'^'^'P) + 1 ■ 

Here Z = e^^. Writing the volume of the momentum element (Pp = Anp'^dp and changing 
of variables [5p^/{2M) — X, these two equations can be written in the forms. 



^^ifJ^flSrM-^.^) (5) 



and 



U = %m.f"^^ J dhU,,{z^^) + ^(^)''7 d\V{r,z)h,,{z.). (6) 
Here = z~^e^'^^'^'^\ and the Fermi integral fi{zr) is given by 



1 rco ^dx 

where r(Z) = /q°° t^~^e~^dt is the Gamma function. For low temperatures (or large values), 
we expand fi{zj.) as an asymptotic series- commonly known as Sommerfeld's lemma 10|, 



fl{Zr) 



\nZr- 



V{l + l] 



! + /(/-! 



TT 



6 (In Zr 



+ /(/-l)(/-2)(/-3) 



360 ilnZrY 



+ ... 



(8) 



Writing the volume element / d^r = 2 Jq° dz Jq^^^^ 27Tr^dr±, with the edges^ of the atomic 
cloud r'j_{z) = \Jb'^/ (47^) + (/i — bzz'^)/'y — b/ (27) and zq = fi/b^ [note b and b^ are defined 



^ Effective chemical potential /i — V{r, z) monotonically decreases from center to the edge of the trap. So 
that the edge of the trap is determined by the condition jj. — V{r, z) — 0. 
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in Eq ([T])], we can perform the spatial integration to obtain analytical expressions for the 
number of atoms and the total energy; 



N 



9>^ 



7r|7| 



OTT 



(9) 



and 



E 



(10) 



171 Vl7l 

where we have defined the dimensionless variables E = U/{huJr), aspect ratio^ of the 
harmonic trap A = ujr/<^z, 7 = ft-l/{M^^r)^ = /^/(^^^r), T = ksT/^hUr), and 
A = (l/16)[l/(|7|/i)]. The dimensionless functions I]\fi{x), I]\f2{x), /_bi(x) = 4J[/i(x)/5 + 
4[/(73(x) + Iu4,{x)]/3, and Ie2{x) = Iu2{x)/2 + [Iu5{x) + Iu6{x)]/6 are defined below. The 
upper sign is for positive values of 7 while the lower sign is for negative values of 7. 



/^i(x) = — ( ± 8(±1 + x^/^TT T v^(15 ± 20a; + 8x'^)7i) 
80 V / 



(11) 



^ , , , ^/x7^ Vil+XTT 
In2[X) = ±— — ± 



(12) 



Iui{x) 



— ( ^ — v^TT (Sx^ ± 16x + 11) T —Vx f±256x^ + 616x^ ± 560x + 175 
96 V 16 ^ ^ 112 ^ 



IT 



48 ( ± 4x^/2 ^ 5^,5/2 ^ 4^3/2 ^ ^ 



7a/xVx±T 



(13) 



1 /2 ^ 2x^/2 + 1 



(14) 



Jj;3(x) = — f 24x^/^ - 24 ± Ix^ ± 56x^/^ ^ 44 ± Ix^ + 35x^/^ 
280 V 



-16Vx±lx±4Va^±l TT 



(15) 



^ The aspect ratio of a thermal cloud is the ratio of the radial to axial size in a harmonic potential. 
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I^Jx) = — (six ± iy/\ - -y/x f 16x^ ± 56x2 ^ ^ \ 
560 V 9. \ 'I 



'C/5(X) = ^ ^ (17) 

6va; ± 1 



2 ^ 3 Jx(x ± 1) ± 4x - 2^^X1x^/2 + 2x2 
lu^ix) = ^ ^ . (18) 

The recurring ± sign originates from the sign of 7 in the integrals. Finally, we derive heat 
capacity C = {dU / dT)\]^ of the system using Eq. ( ITU]) . First, we write 



dU 



_ ke dU 

N huJr df 



N 



-^ + ^-^. (19) 

dT dfidT ^ ^ 



As the number of atoms in the trap is fixed, using the condition dN/dfj, = 0, we find dfi/dT 
to obtain the heat capacity. 



10/iV9^./l^ 167r 7P/2 dA J ImiA) 



III. RESULTS 

In this section we present thermodynamical properties as a function of temperature and 
trap potential parameters. For our calculations, we use a large number of atoms (1.0 x 10^) 
to ensure the validity of our semi-classical treatment. First we solve Eq. [5] for /i for given 
N, T, and A. As can be seen from FIG. [H while the chemical potential increases with the 
anharmonic term, it decreases as the aspect ratio(A) increases. As more energy is needed to 
excite the fermions in low aspect ratio traps, the chemical potential in low aspect ratio traps 
is larger than that of higher aspect ratio traps. As we are using large number of atoms for 
our calculations, anharmonic term cannot be too negative. If the anharmonic term is too 
negative relative to the harmonic confinement, the trapped atoms cannot be confined as the 
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FIG. 1: Chemical potential as a function of anharmonic term 7 for fixed number of atmos N = 
1.0 X 10^. We fixed the temperature to be UbT = huj^- Different lines represent different aspect 
ratios A. 

potential is not bounded from below. The value of the anharmonicity at the point where 
it dominates over the harmonic confinement will be called critical anharmonicity hereafter. 
This critical anharmonicity depends on the number of atoms in the trap. For a set of 
representative parameters, this critical anharmonic value as a function of atom number is 
given in FIG. El Chemical potential for a range of negative anharmonic potential is given 
in FIG. El 

In FIG. m we plot the total energy as a function of anharmonic term. For given values 
of atomic number iV, temperature T, and aspect ration A, we first calculate the chemical 
potential from Eq. IH and then we use this chemical potential to calculate the energy from 
Eq. [ini As can be seen in FIG. HI for fixed number of atoms and fixed aspect ratio A, 
dimensionless energy increases as we increase the anharmonic term 7. This is because of the 
trapping potential energy increases with increasing 7. 

The low temperature behavior of the heat capacity calculated from Eq. (120|) is shown in 
FIGS. [MZl The temperature dependence of the heat capacity for different values of 7 and A is 
shown in FIG. [Hand FIG.O As we have considered up to the linear order in temperature, the 
specific heat is linear with temperature, however, the slope change significantly as we change 
the values of 7 or A. The heat capacity of the polarized Fermi gas increases monotonically as 
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FIG. 2: The critical negative anharmonic dependence on number of atoms in the trap for ksT = hojr 
and A =1. The atomic cloud is unstable below this critical value. 
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FIG. 3: Chemical potential for negative anharmonic potentials. We fixed the temperature, number 
of atoms, and aspect ratio to be ksT = hoJr, N = 1.0 x 10^, and A = 0.5. For 7 < —0.001, the 
system is unstable. 
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FIG. 4: Total energy as function of anharmonic term 7 for fixed number of atmos N = 1.0 x 10^. 
We fixed the temperature to be ksT = hojr and the aspect ratio A = 0.5. For these selected 
parameters, absolute value of critical I7I is 0.001. 



the temperature increases. This behavior is completely different from a Bose gas HI] where 
the heat capacity shows a maximum at supefiuid- normal phase transition point. As can be 
seen from the figures, the heat capacity decreases with increasing anharmonic term 7. This 
decrease is due to the confinement effects which can be clearly seen from Eq. | 



IV. CONCLUSIONS 



We have considered a spin-polarized Fermi gas trapped in an anharmonic potential. 
Applying a semiclassical treatment to the particle number and energy, and using a low- 
temperature series expansion method, we have derived analytical expressions for several 
thermodynamic quantities. We have discussed the trap parameter dependence on these 
quantities and the results were presented as a function of anharmonic trap parameter and 
temperature. The results of our paper can be tested in a lab by a similar experimental 
procedure carried out by Bretin et al. [5i] for Bose-Einstein condensates in a quartic trap. 

In the present work, we considered fully polarized fermions. At low temperatures and 
low densities, s-wave short-range interactions are prohibited and higher wave scattering are 
strongly suppressed. The most natural extension of this work would be the inclusion of 
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FIG. 5: Specific heat as function of temperature for fixed number of atmos N = 1.0 x 10^. We 
fixed the aspect ratio to be A = 0.5. 




f 

FIG. 6: Specific heat as function of temperature for fixed number of atmos N = 1.0 x 10^ and 7 
= 0.5. 

interaction. This can be done by adding another Zeeman state of atoms and consider the 
system as a two component fermions. A particular interesting regime is so called uintarity 
regime, where the scattering length between atoms is infinite. This regime is located near 
the Feshbach resonance 13]. At this point, the interaction parameter disappears from the 
physical quantities. Therefore, all thermodynamic quantities are universal function of the 
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FIG. 7: Specific heat as function of 7 for negative 7's. We fixed the number of atmos iV = 1 x 10^. 
We fixed the ksT = fioJr- 



dimensionless temperature kBT/ep, where ep is the Fermi energy 13|. However, as there is 
no small parameter, it is difficult to reliably calculate thermodynamic quantities close the 
the unitarity regime. Therefore, one has to rely on good experimental data of energy, heat 
capacity, etc. 1^. Then by comparing experimental data at unitarity regime with non- 
interacting calculation of thermodynamic quantities allows one to find the scaling factors. 
An accurate determination of these scaling factors of unitary Fermi gases is one of the 
challenging and fascinating theoretical problem at the moment 15 1. 
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